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We analyze theoretically correlation effects on the transport properties of a benzene molecular 
transistor, that are mediated by interactions between the motion of the nuclei and the transmitted 
charge. We focus on a single phononic mode which allows us to derive an analytic expression for 
the current. The results provide transparent interpretations of various features of the highly non- 
linear current-voltage characteristics, which is experimentally accessible through resonant inelastic 
electron tunneling spectroscopy. The theoretical approach we present is generally applicable to cor- 
related nano-scale devices and considers exactly electron-electron interactions and the coupling to 
vibrational degrees of freedom in the device while the tunneling coupling to the leads is treated 
perturbatively. 



Vibrations of atoms in molecules or solids play an im- 
portant role for many intriguing phenomena in various 
fields of physics. In condensed matter physics novel states 
often emerge from phonon mediated collective behav- 
ior. Among them are conventional superconductivity [1], 
colossal magnetoresistance [2], and thermoelectricity [3]. 
Localized vibrational modes are also believed to stabi- 
lize the a-helix in protein and to be responsible for stor- 
age and transport processes [4]. In molecular electronics, 
the electron-phonon interactions strongly influence the 
conduction properties, contribute significantly to heat- 
ing, and lead to nonlinear effects, including bistability, 
negative differential conductance, and hysteretic behav- 
ior [5-16] . In these systems signatures of the interaction 
between charge carriers and molecular vibrations can be 
faithfully characterized by (resonant) inelastic electron 
tunneling [(R)IET] spectroscopy [15, 17, 18]. The im- 
portance of phonon modes on the conduction properties 
of molecular transistors has been demonstrated in various 
seminal experiments [19-26]. 

A wide range of theoretical studies concerning noncqui- 
librium properties of nanoscale devices are based on the 
nonequilibrium Green's function (NEGF) method eval- 
uated in the frame of the local density approximation 
(LDA) (see e.g. [27, 28] and references therein). On 
the one hand, the LDA-NEGF approach is ab-inito and 
as such allows to study realistic devices, on the other 
hand, correlation effects, irrespective whether they are 
of electron-electron or electron-phonon character, are in- 
cluded only to low order. This can lead to unphysical 
results provided the interaction is strong [29]. 

In the theoretical approach, presented here, we start 
out from the exact expression for the current within the 
NEGF framework and employ a strong coupling cluster 
approximation to evaluated the required Green's func- 
tions (GF) [30, 31]. This allows to exactly take into ac- 
count all interaction of the nanoscale device itself, be 
it of electron-electron or electron-phonon nature, while 
treating the influence of the leads perturbatively. As an 



application, we investigate the nonequilibrium transport 
across a benzene molecular transistor [Fig. 1 (a)] and ob- 
serve that at low temperatures phonon mediated inter- 
actions leave a phononic fingerprint qualitatively similar 
to those observed in (R)IET spectroscopy [19-26]. 

At low temperatures we obtain a stepwise increase 
of the current at voltage differences given by twice the 
phonon energy, Fig. 1 (b). The stepwise increase of the 
current manifests in the conductance and in the differ- 
ential shot noise through multiple, equally-spaced peaks, 
Fig. 1 (c) and (d), as well as in a branching of the dia- 
mond structure in the stability diagram, Fig. 3. At low 
bias voltage the current is significantly suppressed due 
to electron-phonon interactions, which is reminiscent of 
the Franck-Condon blockade in molecular devices [6] , see 
inset of Fig. 1 (b) . 

Theoretical description. — We consider a molecular 
transistor consisting of a single benzene aromatic ring 
coupled to metallic leads, Fig. 1 (a). The benzene 
molecule is modeled by the Hamiltonian 

Hb = H e + H p + H ep . (1) 

The first term describes the motion of the electrons, 

H e = [ _t (4r c i+l<r + h - C + e 4r C «r > 

where t is the hopping amplitude and e the local energy, 
which can be controlled by the gate voltage V g , i.e., e — > 
e + eV g . Electrons with spin a on orbital i are created 
(destroyed) by c\ a (c ia ). The second term H p is the bare 
phonon part, which we transform into the eigenmodes 

m 

where b^ (b m ) create (destroy) phonons in mode m with 
energy huj m . The last term in (1) incorporates the Su- 
Schrieffer-Heeger [32] electron-phonon interaction 

H ep = ]T 9m D^ +1 \bl + b m )(clc J+la + h.c.) 

771 J iff 
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FIG. 1. (Color online) (a) Illustration of the phonon breath- 
ing mode in a benzene molecular transistor (red spheres) 
which is coupled to metallic leads (blue spheres) . (b) The bias 
voltage Vb induces a current j that is significantly modified by 
the interaction between the nuclear motion and the charge. 
At low Vb the current is suppressed [6], see inset, whereas 
at large Vb strong nonlinear behavior leads to a stepwise in- 
crease of the current, also manifesting in (c) the conductance 
and (d) the differential shot noise through multiple equally 
spaced peaks. The data is evaluated at T = 10K. 



with g m := g^J h/2Muj mi g the electron-phonon interac- 
tion, M the mass of the nuclei, and D m the displacement 
matrix of eigenmode m, which is given by 



with 



G {0,1,... 5} 



The metallic leads are described by semi-infinite tight- 
binding chains with hopping t ai on-site energy e Q , and 
chemical potential [i a (a = /eft or right). We assume 
that the bias voltage eV& = jJLi — fi r is applied symmet- 
rically, i.e., fj,i = ei = —fi r = —e r and that the leads are 
equilibrated, i.e., that the occupation follows the Fermi- 
Dirac statistics f a = [e( e- ' la )/ fcBT + l] -1 , where ks is the 
Boltzmann constant and T the temperature. 

In this work we use the NEGF framework to evalu- 



ate the nonequilibrium steady-state properties along the 
lines presented in [30] . Due to the electron-phonon inter- 
action it is in general not possible to solve the problem 
exactly and thus we employ Cluster Perturbation The- 
ory (CPT) [33]: In short, the system is first divided into 
several clusters, each of which can be solved exactly ei- 
ther analytically or numerically, taking correlations on 
the cluster level into account exactly. Strong coupling 
CPT is then employed to connect the individual clusters. 

Here, we divide the molecular transistor into three 
clusters: (i) the central region consisting of the ben- 
zene ring with vibrational degrees of freedom and (ii) 
the left and (iii) right metallic leads which are connected 
to the molecular device. From the exact solution of the 
isolated interacting central region we obtain an approxi- 
mate self-energy of the full system. The results can also 
be systematically improved by incorporating lead sites 
into the central system [30, 31]. Based on considera- 
tions similar to those of Meir and Wingreen [34] we have 
shown [35] that within this CPT approach the current 
is of the Landauer-Biittiker form [36, 37] even though 
correlations are included exactly on the cluster level. In 
particular, we find 



ds 
2^ 



(fi - fr) Tr[T] 



where T is the transmission coefficient matrix 



T:=g R (e)r r (e)g A (e)T l (s) 



(2) 



(3) 



This matrix contains the retarded (advanced) GF Q R 
(Q A = (£«)t) which within CPT is given by (G 11 )- 1 = 
{9 R ) 1 — (Ej + £ r ), where g R is the exact retarded GF 
of the isolated benzene molecule, £ Q :— T ca g R a T ac takes 
into account the effect of lead a on the GF, T CQ is the 
tunneling coupling between the central device and lead a, 
and g R a stands for the retarded lead GF [30]. In addition, 
the transmission T depends on the imaginary part of the 
lead induced self-energy via T a = — 2Im£ a . Similarly, 
the expression for the shot noise, i.e., current-current cor- 
relations, in the zero frequency limit reads [38] 



S(u -> 0) := S 



ds 
2^ 



[/<(l-//) + /r(l-/r)]Tr[7l 



+ (fi-fr) 2 Tr[(l-T)T] 



(4) 



In principle the nonequilibrium properties can be de- 
termined by solving the eigenvalue problem of nu- 
merically taking all phonon modes into account using the 
formulas discussed above. Within this approach it is also 
straight forward to introduce electron-electron interac- 
tions giving rise to Kondo physics provided the molecule 
is (asymmetrically) electron or hole doped [39]. 

In the present work, however, we want to focus on the 
effects generated by the electron-phonon interaction. In 
order to gain insight into the problem we focus on a sin- 
gle phonon mode, in particular on the breathing mode 
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FIG. 2. (Color online) Temperature dependence of (a) the current, (b) the conductance, and (c) the differential shot noise 
versus bias voltage Vj, curves. At temperatures T much lower than the phonon energy ui^ — 0.086eV, which is the smallest 
energy scale in Eq. (1), the phonon mediated steps and peaks, respectively, are pronounced and clearly visible. For increasing 
temperature these features are washed out and at room temperature T = 300K hardly any remain. 



illustrated in Fig. 1 (a), which is the one with the lowest 
excitation energy w M of the isolated benzene ring. An- 
other appealing feature of this mode is the fact that the 
resulting interacting model of the central region can be 
solved analytically and allows to obtain a transparent in- 
sight into the non-equilibrium physics of such a molecular 
electronic device. 

Analytic solution. — For a single phonon mode fi (not 
necessarily the breathing mode) we can readily eliminate 
the electron-phonon interaction by a unitary transfor- 
mation of the Lang-Firsov type [28]. To this end we 
first diagonalize the displacement matrix D^X — XA, A 
is the diagonal matrix with the eigenvalues and X the 
corresponding matrix of eigenvectors, and transform the 
electron operators accordingly d = Ac, where we intro- 



duce the vector notation c 



t = (c\. rt 



-Ll>- 



The 



Lang-Firsov like unitary transformation is generated by 
the operator 



(5) 



All operators that are subject to the unitary transforma- 



tion are denoted by O := e Oe 
results in 

fctr ka 



which in particular 



(6) 



The Lang-Firsov transformation eliminates the explicit 
electron phonon interaction and results in an effective 
Hamiltonian Ht — H e + hu^wb. The electronic GF 
entering the transport properties of the molecular junc- 
tion [Eq. (3)] is transformed accordingly, i.e., g{k,uj) — > 
g(k,uj) by replacing d\ a — > d\ a in its definition. For the 
particular case of the breathing mode the matrices 
and the matrix describing the nearest neighbor hopping 
commute and thus can be diagonalized by a common set 
of eigenvectors. In this case we find = ^/2/3cos k and 
X n k — exp(ikn)/\/6. The d k thus create electrons with 



a given (angular) momentum k in the molecule. In the 
transformed Hamiltonian the electronic part becomes 



H e — J^(e — 2icosfc)nfc cr — U(S^. n ka cos A;) 2 



(7) 



After this transformation, the electron and phonon de- 
grees of freedom are completely decoupled but instead 
an attractive electron-electron interaction proportional 

to U := g^j- emerges. This interacting Hamiltonian 
can be solved exactly. The eigenvectors are of the form 
\n) \Kf) \K±), i.e., they are tensor products of the oscilla- 
tor eigenvectors \n) and the slater determinants \K a ) := 

II* J 



lkeK a 

7t 



dl |0) composed of the molecular orbital opera- 
tors d\ dressed by bond-distortions that are encoded in 
the Lang-Firsov factor of Eq. (6). For the electron GF 
for a single spin-species we find 



g(k,uj) = 



G(k > k F ) 



6(fc < k F ) 



- (Kn - Eo) hio + (E- n - E ) 



(8) 



with ay. — yU cos k and the energies are 

E = eN a - 2tC - UCl 
E^ n = nhuj^ + (Co ± l)e - 2t{C ± cos k) - U(C ± cos k) 2 

Co = cos q ■ 

a <jeFS 

Ao is the equilibrium number of electrons in the benzene 
molecule and FS denotes the Fermi sea. The physical 
properties of the system are particularly obvious in the 
differential conductance which to order 0(V 2 /T 2 ) and 
0{t! 2 jt\) is given by 

dj 



dV b 4for 



[T(V b /2) + T(-V b /2)] 



(9) 
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Along with Eqs. (3) and (8) it can be deduced that cur- 
rent pathways open up, whenever VJ,/2 is equal to an the 
excitation energy ftw = E^ n + huj^n in the particle sector 
or Hw = —E^ n + huj^n in the hole sector. 

Transport properties. — Next, the impact of phonons on 
the transport properties shall be investigated guided by 
the breathing mode of a benzene molecular transistor. 
We use the parameters reported in [13, 40] for the ben- 
zene molecule and the electron phonon coupling; in par- 
ticular we set t = 2.5eV, e = — 1.5eV, hu^ = y/K/M — 
0.086eV, and g^ = 0.18eV. In addition, for the metallic 
leads we choose t a = 12eV close to the wideband limit 
and the coupling between the leads and the benzene ring 
is t! = 0.4eV. Therefore the smallest energy scale in the 
model is the phonon energy Huj^. The current evalu- 
ated at T — 10K (the temperature enters through the 
Fermi-Dirac distribution functions of the leads) for the 
rigid molecule as well as the molecule with the breathing 
phonon mode is shown in Fig. 1 (b). At low bias voltage 
we observe a suppression of the current due to electron- 
phonon interactions, reminiscent of the Franck-Condon 
blockade [6], see inset. Due to the electron-phonon me- 
diated correlations (second term in Eq. (7)) the gap in 
the density of states of the benzene molecule widens con- 
siderably. Consequently, the threshold for the current 
increases to Vb ~ 4eV. Above the threshold the current 
increases stepwise, at voltage differences given by twice 
the phonon energy hui^, until the next plateau is reached. 
The new current pathways are related to the phonon ex- 
citations and the factor 2 is due to the argument 14/2 in 
Eq. (9). 

This also manifests in the differential conductance 
Fig. 1 (c) and in the differential shot-noise (d) through 
pronounced, equally spaced peaks. At the maxima of 
the conductance, the differential shot-noise, correspond- 
ing to the derivative of the current-current correlation 
for uj — > 0, exhibits dips [38]. At low temperatures this 
is observable for a few low-energy peaks of the phonon 
branching. A detailed analysis of the differential con- 
ductance shows that the peaks are near the poles of the 
cluster GF g of the decoupled device. This is reasonable, 
since the effective coupling V/t' 2 is small. Therefore, the 
impact of S Q stems primarily from its imaginary part 
which determines the width and the weight of the indi- 
vidual peaks. So the measurement of these quantities can 
be used as a diagonistic tool to analyze the device and ex- 
tract information about the electronic and the phononic 
excitations. In particular the equidistant peaks in the 
conductance and the tips in the differential shot-noise are 
a fingerprint of the bare phonon excitations. One should 
notice, however, that for general phonon modes this is 
not entirely correct, since the electronic excitations will 
depend on the phonon occupation n, due to H e . 

For increasing ambient temperature the phonon fea- 
tures are smeared out, since the phonon energy lj^ — 
86meV is rather small and comparable to room temper- 
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FIG. 3. (Color online) (a) Current as a function of the bias 
voltage Vb and the gate voltage V g and (b) corresponding sta- 
bility diagram (i.e., differential conductance). The electron- 
phonon interaction increases the area of the central diamond 
and induces a splitting of the lines in the stability diagram, 
(c) Fowler-Nordheim plots for a few selected values of the 
gate voltage V g specified in the legend. Data evaluated for 
T = 10K. 



ature T = 300K = 26meV. In Fig. 2, the current, the 
differential conductance, and the differential shot noise 
are shown at different temperatures T = {10, 77, 300}K. 
With increasing temperature the steps in the current are 
smoothed out and the peaks in the current and the differ- 
ential shot noise decrease significantly. In addition, the 
dips in the shot noise vanish. 

Finally, we study the consequence of orbital gating, 
where the local on-site energy of the benzene molecule 
is changed by a gate voltage, which we restrict to the 
HOMO-LUMO gap of the noninteracting molecule. The 
current as a function of the gate voltage V g and the bias 
voltage Vb exhibits the celebrated diamond structure cen- 
tered around particle-hole symmetric point V g = 1.5V, 
Fig. 3 (a) . As compared to the rigid molecule the area 
of the center diamond is increased due to the electron- 
phonon interactions and the transition between the dia- 
monds involve several steps, best visible in the stability 
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diagram, corresponding to the conductance in the V g -Vb 
plane, shown in Fig. 3 (b) . The Fowler-Nordheim plots 
(c) which show ln(j /V b 2 ) as a function of Vj, -1 compare 
well with the experimental observation in Ref . [26] . From 
these plots different tunneling regimes become apparent: 
At low bias voltage, the "logarithmic" dependence indi- 
cates direct tunneling, whereas for increasing bias volt- 
age the curve has a large negative slope, indicating field 
emission [26]. Close to this bias voltage phonon medi- 
ated interactions manifest as small wiggles in the curves. 
When further increasing the voltage the next plateau is 
reached which is accompanied by a reentrance in the di- 
rect tunneling regime. Similarly to [26] the minimum of 
the curves move to larger bias voltage for increasing gate 
voltage reflecting the diamond shape of Fig. 3 (a) and (b) . 

Conclusions and Outlook. — In summary, we explored 
phonon mediated correlation effects on the transport 
properties of a single benzene-ring transistor. We fo- 
cused on an analytically tractable single-mode model 
that reveals a deeper understanding of the impact of 
phonons on experimental (R)IET spectra. The theo- 
retical approach we presented allows to determine non- 
equilibrium properties of strongly correlated nano-scale 
devices. When solving the interacting device totally nu- 
merically, it also allows to include other phonon modes 
and electron-electron interactions. 
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